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I. INTRODUCTION 



Recent cosmological and astrophysical data gathered 
from the observations of supernovae SNe la \E, cosmic 
microwave background radiations via WMAP Q, galaxy 
redshift surveys via SDSS [3| and galactic X-ray |4| con- 
vincingly suggest that the observable Universe experi- 
ences an accelerated expansion phase. It is well-known 
that the simplest and elegant way to explain this behav- 
ior is the inclusion of Einstein's cosmological constant 
[5(, however the two deep theoretical problems (namely 
the "fine-tuning" and the "coincidence" one) led to the 
notion of 'dark energy'. The dynamical nature (i.e. com- 
position and origin) of dark energy, at least in an ef- 
fective level, can arise from various scalar fields, such 
as a canonical scalar field (quintessence) [6|, a phantom 
field Q, that is a scalar field with a negative sign of 
the kinetic term, or the combination of quintessence and 
phantom in a unified model named quintom [8(. One of 
the first quintom works which appeared is [9| . Recent re- 
view on dynamical DE from modified gravity perspective 
can be found in [10] . Also, there are other momentous 
works on dynamics of the scalar models for accelerating 
expansion phase of the universe with different methods 
and terminologies (llj. Recently, a new type of the scalar 
models have been investigated more than others [l2j . 

One of the long-standing problems in the standard Big 
Bang cosmology is the initial singularity from which all 
matter and energy originated. Standard cosmology based 
on general relativity offers no resolution to this problem, 
however a recent quantum gravitational model of loop 
quantum gravity (LQG) offers a nice solution. The the- 
ory and principles of LQG when applied in the cosmo- 
logical framework creates a new theoretical framework 
of loop quantum cosmology (LQC). The effects of loop 
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quantum gravity can be described in two possible ways: 
the first one is based on the modification of the behavior 
of the inverse scale factor operator below a critical scale 
factor. This approach has been used to study quantum 
bounces, avoidance of singularities and to produce infla- 
tionary expansion 20]. A second approach is to add a 
term quadratic in density to the Friedmann equation. In 
LQC, the non-perturbative effects lead to correction term 
—p 2 / Pc to the standard Friedmann equation. With the 
inclusion of this term, the Universe bounces quantum me- 
chanically as the energy density of matter-energy reaches 
the level of p c (order of the Planck density). Thus the 
LQC is non-singular by producing a bounce before the 
occurrence of any potential singularity and hence transi- 
tions from a pre-Bang and after-Bang are all well-defined. 
The observational constraints due to the quadratic term 
in (1) are discussed in the literatures [21( where it is 
shown that the model with quadratic correction to den- 
sity is consistent with the observational tests. Thus we 
should not worry on the solar system tests in this model. 
For the Universe with a large scale factor, the first type of 
modification (first approach) to the effective Friedmann 
equation can be neglected and only the second type of 
modification (second approach) is important [20j. Thus 
the dynamics of dark energy have been investigated re- 
cently in LQC using second approach [22 ] . 

In this paper, we address the problem of cosmic co- 
incidence problem in a modern theoretical framework of 
loop quantum cosmology. Here we assume a non-minimal 
coupling between dark energy, dark matter and the un- 
particle component. Since the nature of dark energy and 
dark matter is still unknown, it is possible to have non- 
gravitational interactions between them. We are unsure 
of the form of the interaction, hence we introduce the in- 
teraction terms only phenomenologically. We construct 
a system of dynamical equations containing three equa- 
tions for the three components. We convert them to 
dimensionless form and perform stability analysis. We 
construct four toy models and show that the dynamical 
equations have two possible attractor solutions i.e. dark 
matter dominated and dark energy dominated. Other 



dynamical systems are unstable i.e. one model predicts 
that everything decays and the Universe gets emptier 
void of energy. Another model predicts that the Uni- 
verse will contain only dark energy. 

The model in which dark energy interacts with two 
different fluids has been investigated in the literature. In 
[13j . the two fluids were dark matter and another was 
unspecified. However, in another investigation [14| |. the 
third component was taken as radiation to address the 
cosmic-triple-coincidence problem and study the gener- 
alized second law of thermodynamics. In this article, we 
choose the third component as unparticle (l5| . following 
}16| thereby generalizing their study from the general rel- 
ativistic cosmology to the loop quantum cosmology. An 
unparticle is based on the hypothesis that there could 
be exact scale invariant hidden sector resisted at a high 
energy scale. The fundamental energy scale of unparticle 
is far beyond the reach of today's accelerators, there is 
a possibility that this new unparticle sector could affect 
the low-energy phenomenology. An interesting feature of 
unparticle is that it does not have a definite mass and in- 
stead has a continuous spectral density as a result of scale 
invariance. Moreover, the equation of state of unparticle 
w u is positive unlike dark energy and it interacts weakly 
with standard model particles. We consider the question 
how the evolution of Universe is affected when the unpar- 
ticle takes part in the interaction with dark energy and 
dark matter. 

This paper is organized as follows: in section II, we 
construct a cosmological scenario in which dark energy 
interact with dark matter and unparticle. Then we 
present the formalism of autonomous dynamical system 
which is suitable for phase space stability analysis. In 
section III, we study the stability of the dynamical sys- 
tem by choosing different coupling functions. Finally we 
briefly discuss our results in the last section. 



II. OUR MODEL 

Applying the techniques of loop quantum gravity to 
homologous and isotropic spacetime leads to the so-called 
loop quantum cosmology. Due to quantum corrections, 
the Friedmann equations get modified. The big bang 
singularity is resolved and replaced by a quantum bounce 
[17j . For a brief summary on loop quantum cosmology, 
see [18j. Considering quantum correction, the modified 
Friedmann equation turns out to be (in the case of k = 0) 
181 
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where p = p m + pd + pu, where p m , p d and p u represent 
the energy densities of matter, dark energy and the radi- 
ation. Also p c = .„ 2 3 /9pi determines the loop quantum 
effects, ppi is the Planck density and 7 is the dimension- 
less Barbero— Immirzi parameter. This parameter could 



be fixed as 0.2375 in order to give the area formula of 
black hole entropy in loop quantum gravity [l9| . The ob- 
servational constraints due to the quadratic term in (1) 
are discussed in the literature [2l| , where it is shown that 
the model with quadratic correction to density is consis- 
tent with the observational tests. Thus we should not 
worry on the solar system tests in this model. 

Due to the corrected term in ([1]), the big bang sin- 
gularity is replaced by a quantum bounce happening 
at p c . The bounce is supposed to happen when the 
matter— energy density reaches the critical value p c . 
However the numerical simulations show that modified 
Friedmann equations are valid in the whole cosmic evo- 
lution including the bounce |17| . 

Another FRW equation is 
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For a spatially flat Universe, the total energy conserva- 
tion equation is 



p + 3H(p + p)=Q, 



(3) 



where H is the Hubble parameter, p is the total energy 
density and p is the total pressure of the background 
fluid. 

We assume a three component fluid containing matter, 
dark energy and unparticle having an interaction. The 
corresponding continuity equations are [161 ] 



pd + 3H(p d +p d ) = Ti, 

pm + 3Hp m = r 2 , 

p u + 3H(p u +p u ) = r 3 , 



(4) 



which satisfy collectively (|3|) such that Ti + T 2 + T 3 = 0. 
We define dimensionless density parameters via 
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Making use of parameters in ([5]) in the modified Fried- 
mann equation ([1} is 



(x + y + z)(l 
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Using |T]) and ©, we can write 



#3 / w d x + w u z\ 

H z 2 \ x+y+z I 

The equation of state parameter of the total fluid is 

p w d x + w u z 

Wtot = - = : : ■ 8 

p x+y+z 

The continuity equations Q in dimensionless variables 



reduce to 



III. ANALYSIS OF STABILITY IN PHASE 
SPACE 
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Above the primes denote differentiation with respect to 
N = In a. The coupling functions r,, i = 1,2,3 are in 
general functions of the energy densities and the Hubble 
parameter i.e. Ti(Hpi). The system of equations in (|9|) 
is analyzed by first equating them to zero to obtain the 
critical points. Next we perturb ([9]) up to first order 
about the critical points and check their stability. 

Following [161 ] , one can consider various coupling func- 
tions to model the interaction. However all the models 
described here are phenomenological. There is, at yet, 
no deep theoretical justification for assuming this inter- 
action in the absence of the quantum gravity. There has 
been a striking resurgence in modeling dark energy via 
such interacting models as they explain the astrophys- 
ical data with the desired accuracy. Moreover, in such 
models, one can obtain attractor solutions (stable equi- 
librium points against perturbation) at which the energy 
densities come at equilibrium and thereby explain cosmic 
coincidence. These models have very little fine tuning 
problems since these contain only one arbitrary coupling 
parameter. In literature, this parameter has been con- 
strained via astrophysical data from supernovae, cosmic 
background radiation, baryon acoustic oscillations, gas 
mass fraction in galaxy clusters, the history of the Hub- 
ble function, and the growth function [231 ] . The signature 
of the coupling parameter is of central importance: its 
positive (negative) sign gives the direction of interaction. 

There are some drawbacks as well for the interacting 
dark energy models: they are not clearly distinguishable 
from the non-interacting ones, in the light of the obser- 
vational data [24]]. Another drawback is that there is 
arbitrariness in the choice of the interaction term (i.e. 
products of Hubble parameter with the energy densities) ; 
there are models in which only the coupling parameter or 
the dark energy state parameter is taken as a function of 
scale factor to model interaction without employing the 
densities or the Hubble parameter [25j . 



In this section, we will construct four models by choos- 
ing different coupling forms Ti and analyze the stability 
of the corresponding dynamical systems about the crit- 
ical points. We shall plot the phase and evolutionary 
diagrams accordingly. 



A. Interacting Model I 

We consider the model with the following interaction 
terms 



Ti = —QbHp x , T2 



3bHp x 



(10) 



where 6 is a coupling parameter and we assume it to be a 
positive real number of order unity. Thus (ITU1) says that 
both matter and unparticles have increase in energy den- 
sity with time, while dark energy loses its energy density. 
Therefore, it is a decay of dark energy into matter and 
unparticle. 

Using (fTDj). the system ([9]) takes the form 
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w d x + w u z 
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x + y + z 
—3x — 3wdX — 66a;, 

, / WdX + WuZ\, 
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x + y + z 

—3y + 3bx, 

I o (1 , w d x + w u z 
z = 3z(H 

V x + y + z 

— 3z — 3w u z + 3bx. 
There are three critical points: 

• Point A x : (0,0,1), 

• Point Si : (0,1,0), 



. Point d = ( {w f 2 % w h d+2 \T-l 

1 V w d (w d -\-2b—w u ) -\-w u b ' 



(ii) 
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b(w d +2b-wy) 
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The eigenvalues of the Jacobian matrix for these critical 
points are: 



• Point A\ : X\ — w u — Wd — 26, A2 
A3 = -1 -w u , 
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• Point Bi : Xi = -(w d + 26), A 2 = -1, A 3 = -w u , 

• Point Ci : Ai = -(1 + w d + 26), 

X _ _,„ ( (w d +2b)(w d +2b-w u ) \ 
A2 w d,y Wd ( Wd+2b - Wu ) +Wub J, 

\ _ { (w d +b)(w d -w u ) A 
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FIG. 1: Variety of x c , y c and z c with b at the critical point 
C'i for fixed Wd = —1.2 and w u — 0.28. Here the coupling 
constant b is located in the region Wd < —b. 



FIG. 2: Variety of x,y,z as a function of the N = ln(a). 
The initial conditions chosen are x(0) = y(0) — z(0) = 0.3, 
w u — 0.28 and Wd = —1.2 and b — 0.5. We observe that the 
dominant part of in Model I is matter. 



For the point Ai , the eigenvalue A2 is non- negative, which 
indicates that A\ is not a stable point. For the point B±, 
we also find that both of the eigenvalues Ai and A3 are 
negative when Wd > —26. Therefore the point B\ is the 
stable point. For the point C\ , when Wd < —b and Wd < 
w u , all of the eigenvalues Ai, A2 and A32 are negative, 
which means that C\ is a stable point. As figure 1 shows 
below, the stable critical point C\ gives us the intuitive 
picture where the x c (representing DE) decreases while 
y c (representing matter) and z c (representing unparticle) 
increases. Moreover, the situation arises where matter 
dominates asymptotically and unparticle density comes 
after matter. Since C\ is an attractor solution, it implies 
that the dynamical equations ([9]) yield this behavior for 
generic initial conditions. It is interesting to note that 
similar situation arises in Einstein's gravity as well |16| . 
The stable region of C\ is not affected by the EoS of the 
unparticle, but the position of C\ in the phase space is 
decided together by Wd, w u and the coupling constant b. 
From <JSj> , we also learn that the effective total EoS at 
point C\ is 



Wtot 



w d (2 b(w d + 2b- w u ) + w d (w d + 2 b - w u )) 
Wd(w d + 2b- w u ) + bw u 
w u b(2b + w d ) 



Wd{w d + 26- w u ) - bw u 
2 b{w d + 26- w u ) + w d (w d + 26- w u ) 
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The critical point C\ denotes that the DE, DM and un- 
particle can be coexisted in the late-times of the Universe. 
Figure (2) shows the evolution of the functions (x, y, z) 
for a variety set of parameters . As we observe, for a 



set of density parameters Wd = —1-2, w u = 0.28, 6 = 0.5, 
matter is the dominant portion and the unparticle's den- 
sity always remain below the dark energy. One interest- 
ing feature is that for large values of scale factor (the 
smaller values of the rcdshift z) the two different parts of 
the matters in the Universe (unparticle and dark energy) 
tend to the same value and after reaching to this point, 
the dominant part is the matter field. 

Now we change the parameters to a new set Wd = —1-7, 
w u = 0.35, 6 = 0.25. The behavior of the functions (den- 
sities) are very different. Now, if the evolution begins 
from a very large negative red shift, the dominant part of 
the model for all times is the dark energy. As the previous 
case, the unparticle's portion remains under dark energy 
and matter. For some values of large scale factor, the 
three different parts of the matters fields (dark energy, 
matter, unparticle) reach to the same asymptotic's value. 
These evolutionary scheme has been shown in figure (3) 
As we observe in figure (4), for a set of density parame- 
ters — 1 < Wd < —-g,w a > 0, 6 = 0.5, w tot > 0, matter is 
the dominant portion and the unparticle's density always 
remain below the dark energy. It is consistent with the re- 
sult shown in figure (2) in which as the coupling is strong 
enough the Universe will be dominated by DM. Figure (5) 
shows the phase diagram of interacting dark energy with 
DM and unparticle in loop quantum cosmology through 
the coupling terms. The point C\ is the critical point. 
Here we choose the values Wd = —1-7, w u = 0.28, 6 = 0.5 
in the stable region Wd < —b and Wd < w u - 



B. Interacting Model II 

We study another model with the choice of the inter- 
action terms 

r a = -3bH Pd , r 2 - 3bH{p d - Pm ), r 3 = 3bH Pm . (13) 
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FIG. 3: Variety of x, y, z as a function of the N = ln(a). As 
we observe, the dominant part of the model is dark energy. 
The initial conditions chosen are x(0) = y(0) — z(0) — 0.3, 
w u = 0.28 and Wd — —1.2 and b — 0.5. 




FIG. 4: Variety of w tot as a function of the Wd,w u . AS we 
observe, w to t > 0. 



This model effectively describes the situation when dark 
energy loses energy density to matter while the unparticle 
density increases due to interaction with the matter. 
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FIG. 5: The phase diagram of interacting dark energy with 
DM and unparticle in loop quantum cosmology through the 
coupling terms. Here we choose the values Wd = — 1.7, w u = 
0.28, b — 0.5 in the stable region Wd < —b and Wd < w u . 
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The eigenvalues of the Jacobian matrix for these criti- 
cal points are: 

• Point A 2 : Ai = w u — w d — b, X 2 = w u - b, 
A3 = -(1 + w u ), 

• Point B 2 :\i = -w d , A 2 = -(1 + b), 
\ 3 = ±(w u -l-2b) + b, 



• Point C 2 : Ai = -(1 + w d + b), 



\ _ „,, ( w d {b+w d -w u ) 



Aa 



w d —w u )+bw u J ' 
(w d -b)(w d -w u ) 
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For point A 2 , the eigenvalue Aj < if b > w u — w d for 
w d < 0, hence point A 2 is a stable critical point. The 
second critical point B 2 is unstable since Ai > 0. It is 
easy to see that the third critical point C 2 is stable when 
we have the next set of inequalities 



- (1 + u>d) <b < w d 

w d . . 

-w d <b,b< — (w u - Wd) 

w u 



(15) 
(16) 



Figure (3) shows that all the three components start with 
equal energy densities, they evolve differently. Here the 
dark energy dominates over matter and unparticle at 
late times. The energy density of dark energy evolves 




FIG. 6: Variety of x c , y c and z c with b at the critical point 
B2 for fixed Wd = —1.2 and w u — 0.28. 




FIG. 7: Variety of x c , y c and z c with b at the critical point 
Ci for fixed Wd — —1.2 and w u = 0.28. Here the coupling 
constant b is located in the region b < Wd,Wd < w u . 



to x ~ 0.64 and y ~ 0.23 compatible with the obser- 
vations, while the unparticle density falls to zero at an 
early epoch N ~ 2. Clearly this phenomenological model 
of dynamical interaction explains the present state of the 
Universe. The next figure (8), shows, the phase diagram 
of interacting dark energy with DM and unparticle in 
loop quantum cosmology through the coupling terms . 
The point Ci is the stable critical point. Here we choose 
the values Wd = —1-7, w u = 0.28, b = 0.5 in the stable re- 
gion -(1 + Wd) <b <w d ,w u -w d <b,b < %±(w u —Wd). 



C. Interacting Model - III 

The coupling terms containing the product of energy 
densities have been studied previously with the intent to 
reveal some new behavior of the dynamical systems [20j . 




FIG. 8: The phase diagram of interacting dark energy with 
DM and unparticle in loop quantum cosmology through the 
coupling terms. Here we choose the values Wa = — 1.7, iy« = 
0.28, b = 0.5 in the stable region —(1 + Wd) < b < Wd, W u — 
Wd < b, b < ^f-(w u - w d ). 



Let us take the interaction terms [16j 

r a = -Qb^H^papu, r 2 = r 3 = 2bK 2 H- l PdPu . (17) 

The system in ([9]) takes the form 
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There are five critical points: 

• Point A 3 : (1,0,0), 

• Point B 3 : (0,1,0), 

• Point C 3 : (0,0,1), 

• Point D 3 
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Now we must diagonalize the Jacobian matrix near these 
critical points. For points A 3l B 3 ,C 3 we have 
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FIG. 9: Variety of x,y,z as a function of the N = ln(a). 
The initial conditions chosen are ir(0) = j/(0) = 2(0) = 0.3, 
w u = 0.28, w d = -1.2 and b = 0.5. 



FIG. 10: Variety of x,y,z as a function of the N = ln(a). 
The initial conditions chosen are x(0) = y(0) = z(0) — 0.3, 
w u = 0.28, w d = -1.2 and b = 0.5. 



• Point A3 : Ai = — 3(1 + Wd), A2 = 3wd, A3 = 3(36- 

W u +W d ), 



Point B 3 : A x = —3, A2 = — 3w u , A 3 = —'iwj 



Point C 3 : Ai = -3(1 + w u ), A 2 = 3w u ,\ 3 
3(-66 + w u - Wd), 



For point A3 since always w u > 0,Wd < —1 thus this 
point is un stable. Similarly, the point B^^C^ both are 
unstable. The analysis of stability for points D3, E3 are 
so complicated. Indeed , the Jacobian matric in these 
cases, are not diagonal and the behaviors of the eigenval- 
ues are not trivial. Theoretically, we cannot distinguish 
between these points as the stable or un stable points. 
Thus, it is computationally and analytically impossible 
to analyze these cases. Using the same initial condi- 
tions, Figure (9) shows that dark energy density rises 
while matter density dominates over dark energy and un- 
particle till N ~ 1. For N > 1, dark energy density rises 
indefinitely (behaving like phantom energy) while matter 
and unparticle density tends to zero. 



D. Interacting Model - IV 



Consider another model with the interaction terms ll( 
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The system in @ takes the form 

, „ A , w d x + w u z 
x = 6x I 1 H 

V x + y + z 

—3x — iw^x — 9bxz 

WdX + w u z 



y' = 3y(l + 



x + y + z 
3y + 9b(xz — yz), 



(2-x-y-z) 

(2-x-y- z) 



(20) 



Z> = 3z(l+ WdX + WuZ )(2-x-y-z) 
V x + y + z I 

— 3z — 3w u z + 9byz. 
There are six critical points: 

• Point A 4 : (1,0,0), 

• Point B 4 : (0,1,0), 

• Point d : (0,0,1), 

• Point Di : fo, 



3 b 



3ft 



• Point E 4 : U -^+36+^ ) 

1 -2w d w u +w u 2 +3w d b+w d 2 -3wyb 1 w d (w d -w u ) 
3 h(-«>„+36) '3 b(-iu u +3b) J' 



• Points: (i^-1,1^,-^ 



(19) 



Now we must diagonalize the Jacobian matrix near these 
critical points. For points A4, B4, C4 we have 

• Point A4 : Ai = —3(1 + Wd), A 2 = 3wd, 
A3 = 3(w d -w u ), 

• Point B4 : Ai = —3, A 2 = — 3wd, A3 = 3(36 — w u ), 

• Point C4 : Ai = —3(1 + w u ), A 2 = 3(w u — 36), 
A 3 = 3(-36 + ui„ - wo). 
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For point A4 since always w u > 0,Wd < — 1 thus this 
point is un stable, specially for the cross line Wd = — 1 
which in this case, the system is unstable near A±. Sim- 
ilarly, the point -B4 is unstable. But the point C4 is a 
stable point if Wd < — 1, and b > (w u — Wd)/3- The anal- 
yses of stability for points -D4 and E4, are as complicated 
as the Model III. Indeed, the Jacobian matric in these 
cases, are not diagonal and the behaviors of the eigenval- 
ues are not trivial. Theoretically, we cannot distinguish 
between these points as stable or unstable points. Thus, 
it is computationally and analytically impossible to ana- 
lyze these cases. However, the behavior of the dynamical 
equations in (20) is plotted in Figure (10) which shows 
that dark energy density rises rapidly till N ~ 2.5 after 
it decreases sharply (behaving like quintessence), while 
matter and unparticle density decrease and approach to 
zero at N > 2. It is clear that dark energy does not 
decay into matter or unparticle, otherwise their densities 
would have increased. Hence the decay of dark energy 
into some mysterious component is not clear. 



IV. DISCUSSION 

In this paper, we have studied the dynamical behaviors 
when dark energy is coupled to dark matter and the un- 
particle in the background of flat FRW spacetime. The 
analysis was performed in the theoretical framework of 
loop quantum cosmology which is one of the emerging 
fields of quantum cosmology. It was assumed that the 
correction terms due to LQC remain benign and become 



effective near the Planck density. To address the coinci- 
dence problem, we introduced a phenomenological inter- 
action between dark energy, dark matter and unparticle. 
We constructed four toy models: here Model I describes 
a stable attractor solution for a dark matter dominated 
Universe; Model II also contains a stable attractor solu- 
tion for a dark energy dominated Universe; Model III and 
IV are dynamically unstable since either they give com- 
pletely dark energy filled Universe or completely devoid 
of it. 

Comparison of our LQC dynamical models with those 
of Einstein cosmology [16[ reveals the following differ- 
ences: the form of most of the critical points (and hence 
eigenvalues) are different, consequently the conditions 
and regions of stability of these critical points also differs 
from Einstein cosmology. Also a comparison of Fig. 1 
above with Fig. 2 of [16[ shows that the evolution of en- 
ergy densities is much slower in LQC than in Einstein's 
case, however, the asymptotic evolution is the same. Fi- 
nally, unlike [161 ] where all the considered models were 
stable, we have only two stable scaling solutions, namely 
model I and II. In other words, LQC results are very 
different from Einstein's relativistic cosmology. 
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